library(VennDiagram)

## Loading required package: grid

## Loading required package: futile.logger

library(Seurat)

## Loading required package: ggplot2

## Loading required package: cowplot

## 
## Attaching package: 'cowplot'

## The following object is masked from 'package:ggplot2':
## 
##     ggsave

## Loading required package: Matrix

markers <- read.table('../data/markergenes/markergenes-crossed', header=T)

Marker genes of the intersection of Supra and Peri. Supra has 99 marker genes, Peri 71. They have 9 genes in common.

genes.peri <- markers[markers$cluster == 'Peri.white', 'gene']
genes.supra <- markers[markers$cluster == 'Supra.white', 'gene']
peri.supra <- intersect(genes.peri, genes.supra)

grid.newpage()
draw.pairwise.venn(length(genes.peri), length(genes.supra), length(peri.supra), category=c('Peri', 'Supra'), lty=rep('blank', 2), fill=c('#f67770', '#1fbfc3'), scaled=T, cat.fontfamily = rep("helvetica", 2))

## (polygon[GRID.polygon.11], polygon[GRID.polygon.12], polygon[GRID.polygon.13], polygon[GRID.polygon.14], text[GRID.text.15], text[GRID.text.16], text[GRID.text.17], text[GRID.text.18], text[GRID.text.19])

Marker genes of the intersection of Subq and Visce. Subq has 65 marker genes, Visce 22. They have only 2 genes in common.

genes.visce <- markers[markers$cluster == 'Visce.brown', 'gene']
genes.subq <- markers[markers$cluster == 'Subq.brown', 'gene']
visce.subq <- intersect(genes.visce, genes.subq)

grid.newpage()
draw.pairwise.venn(length(genes.visce), length(genes.subq), length(visce.subq), category=c('Visce', 'Subq'), lty=rep('blank', 2), fill=c('#c680fc', '#7dac1f'), scaled=T, cat.fontfamily = rep("helvetica", 2))

## (polygon[GRID.polygon.20], polygon[GRID.polygon.21], polygon[GRID.polygon.22], polygon[GRID.polygon.23], text[GRID.text.24], text[GRID.text.25], text[GRID.text.26], lines[GRID.lines.27], text[GRID.text.28], text[GRID.text.29])

all10x <- readRDS('../data/10x')

Markers intersection Peri and Supra

print(peri.supra)

## [1] "TM4SF1"   "MTRNR2L8" "LY6K"     "MGST3"    "FHL2"     "COL6A2"  
## [7] "TAGLN"    "HSPA8"    "PENK"

markers[which(markers$gene %in% peri.supra),]

VlnPlot(all10x, features.plot=toupper(c('TM4SF1', 'MTRNR2L8', 'LY6K', 'MGST3', 'FHL2', 'COL6A2', 'TAGLN', 'HSPA8', 'PENK')), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)

Markers intersection Subq and Visce

print(visce.subq)

## [1] "GREM1" "HOXB7"

print(markers[which(markers$gene %in% visce.subq),])

##         cluster p_val avg_logFC pct.1 pct.2 p_val_adj  gene
## 190 Visce.brown     0 0.3825937 0.771 0.691         0 GREM1
## 205 Visce.brown     0 0.2671629 0.540 0.324         0 HOXB7
## 220  Subq.brown     0 0.5615335 0.888 0.691         0 GREM1
## 252  Subq.brown     0 0.2726722 0.601 0.324         0 HOXB7

VlnPlot(all10x, features.plot=toupper(c('GREM1', 'HOXB7')), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)

Top markers Peri

print(markers[markers$cluster == 'Peri.white',][1:10,])

##       cluster p_val avg_logFC pct.1 pct.2 p_val_adj     gene
## 1  Peri.white     0 0.8670990 0.928 0.407         0   TM4SF1
## 2  Peri.white     0 0.8098424 0.963 0.681         0     SRGN
## 3  Peri.white     0 0.6577277 0.564 0.161         0     NEFM
## 4  Peri.white     0 0.6504400 0.830 0.468         0     MEST
## 5  Peri.white     0 0.6249721 0.932 0.628         0     RGS4
## 6  Peri.white     0 0.6177285 0.961 0.817         0     PTX3
## 7  Peri.white     0 0.5764050 0.698 0.354         0 MTRNR2L1
## 8  Peri.white     0 0.5387147 0.967 0.802         0   GLIPR1
## 9  Peri.white     0 0.5092244 0.994 0.964         0    CRYAB
## 10 Peri.white     0 0.5058522 0.917 0.656         0     CTSC

VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Peri.white',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)

Top markers Supra

print(markers[markers$cluster == 'Supra.white',][1:10,])

##        cluster p_val avg_logFC pct.1 pct.2 p_val_adj      gene
## 81 Supra.white     0 1.7650025 0.786 0.226         0    IGFBP5
## 82 Supra.white     0 0.9351586 0.977 0.832         0    IGFBP3
## 83 Supra.white     0 0.7779768 0.978 0.800         0    AKAP12
## 84 Supra.white     0 0.6136319 0.804 0.475         0 TNFRSF11B
## 85 Supra.white     0 0.5950695 0.955 0.719         0       DCN
## 86 Supra.white     0 0.5918296 0.884 0.539         0     IFI27
## 87 Supra.white     0 0.5818591 0.808 0.448         0    EFEMP1
## 88 Supra.white     0 0.5744373 0.501 0.176         0      PENK
## 89 Supra.white     0 0.5733726 0.995 0.892         0     MFAP5
## 90 Supra.white     0 0.5585209 0.905 0.656         0    CYP1B1

VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Supra.white',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)

Top markers Visce

print(markers[markers$cluster == 'Visce.brown',][1:10,])

##         cluster p_val avg_logFC pct.1 pct.2 p_val_adj          gene
## 189 Visce.brown     0 0.4954102 1.000 0.999         0          FTH1
## 190 Visce.brown     0 0.3825937 0.771 0.691         0         GREM1
## 191 Visce.brown     0 0.3562843 1.000 1.000         0           FTL
## 192 Visce.brown     0 0.3431171 0.588 0.457         0           LUM
## 193 Visce.brown     0 0.3366891 0.999 0.996         0          GUK1
## 194 Visce.brown     0 0.3158248 0.595 0.430         0 RP11-173B14.5
## 195 Visce.brown     0 0.3156858 1.000 0.998         0          RPS2
## 196 Visce.brown     0 0.3098157 0.930 0.871         0         PLAC9
## 197 Visce.brown     0 0.3089736 0.546 0.365         0         HOXA5
## 198 Visce.brown     0 0.3059777 0.992 0.985         0        IGFBP4

Interestingly, some genes seem to be found because Supra4 and Subq4 are different.

VlnPlot(all10x, features.plot=as.vector(markers[markers$cluster == 'Visce.brown',][1:10,'gene']), group.by='sample_name', point.size.use=-1, nCol=2, x.lab.rot=T, size.x.use=10)

Top markers Subq

print(markers[markers$cluster == 'Subq.brown',][1:10,])

##        cluster p_val avg_logFC pct.1 pct.2 p_val_adj       gene
## 213 Subq.brown     0 1.1149278 0.897 0.722         0      THBS1
## 214 Subq.brown     0 0.9239007 0.974 0.948         0      TIMP3
## 215 Subq.brown     0 0.8931719 1.000 1.000         0      RPS29
## 216 Subq.brown     0 0.8099464 0.402 0.094         0     BCYRN1
## 217 Subq.brown     0 0.7024385 0.831 0.729         0       DKK1
## 218 Subq.brown     0 0.6638062 0.975 0.980         0     RPL36A
## 219 Subq.brown     0 0.6095757 0.973 0.961         0      ATP5I
## 220 Subq.brown     0 0.5615335 0.888 0.691         0      GREM1
## 221 Subq.brown     0 0.5231598 0.654 0.644         0 AC009501.4
## 222 Subq.brown     0 0.5113614 0.997 0.996         0      RPS10

Most of the markers for Subq seem to be found because Subq_3 stands out (RPS29, BCYRN1, RPL36A, ATP5I, AC009501.4, RPS10). Mostly ribosomal genes.

BCYRN1 = brain cytoplasmic RNA